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Abstract 

Extending the parameter space of the three-dimensional (d = 3) Ising model, we search for a 
regime of eliminated corrections to finite-size scaling. For that purpose, we consider a real-space 
renormalization group (RSRG) with respect to a couple of clusters simulated with the transfer- 
matrix (TM) method. Imposing a criterion of "scale invariance," we determine a location of the 
non-trivial RSRG fixed point. Subsequent large-scale TM simulation around the fixed point reveals 
eliminated corrections to finite-size scaling. As anticipated, such an elimination of corrections 
admits systematic finite-size-scaling analysis. We obtained the estimates for the critical indices as 
v = 0.6245(28) and y^ = 2.4709(73). As demonstrated, with the aid of the preliminary RSRG 
survey, the transfer-matrix simulation provides rather reliable information on criticality even for 
d = 3, where the tractable system size is restricted severely. 
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I. INTRODUCTION 



The transfer-matrix method has an advantage over the Monte Carlo method in that 
it provides information free from the statistical (sampling) error and the problem of slow 
relaxation to thermal equilibrium. On one hand, the tractable system size with the transfer- 
matrix method is severely limited, because the transfer-matrix size increases exponentially 
as the system size N enlarges; here, N denotes the number of spins constituting a unit of 
the transfer-matrix slice. Such a limitation becomes even more serious in large dimensions 
(d > 3). Actually, for large d, the system size N(= L d ~ l ) increases rapidly as the linear 
dimension L enlarges, and it soon exceeds the limit of available computer resources. Because 
of this difficulty, the usage of the transfer-matrix method has been restricted mainly within 
d = 2. 

In this paper, we report an attempt to eliminate the finite-size corrections of the d = 3 
Ising model by tuning the interactions parameters. As anticipated, such an elimination of 
corrections admits systematic finite-size-scaling analysis of the numerical data with restricted 
system sizes. To be specific, we consider the d — 3 Ising ferromagnet with the extended 
interactions, 

H = —Jnn SiSj — Jnnn SiSj — Jn SiSjSkSi, (1) 

(id) «M» [idM 

where the Ising spins Si = ±1 are placed at the cubic-lattice points specified by the in- 
dex i] The summations Yltij)> an d Yl\ij,ki] run over a ^ nearest-neighbor pairs, 
next-nearest-neighbor (plaquette diagonal) spins, and round-a-plaquette spins, respectively. 
Within the extended parameter space (Jnn, Jnnn, Jn), we search for a regime of eliminated 
corrections to scaling. For that purpose, we consider a real-space renormalization group for 
a couple of clusters, whose thermodynamics is simulated with the transfer-matrix method; 
see Fig. ^ We then determine a location of the renormalization-group fixed point. Fol- 
lowing this preliminary renormalization-group survey, we perform extensive transfer-matrix 
simulation around this fixed point. Thereby, we show that the corrections-to-scaling be- 



havior is improved around t 



transfer-matrix method 



re fixed point. Here, we utilized an improved version of the 



and succeeded in treating a variety of system 
sizes N — 5, 6, . . . , 15; note that conventionally, the tractable system sizes are restricted to 
AT = 4, 9, 16, ... . Apparently, such an extension of available system sizes provides valuable 
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information on criticality. 

In fairness, it has to be mentioned that our research owes its basic idea to the following 
pioneering, studies: First, an attempt to eliminate the finite-size corrections was reported 
in Ref. where the authors investigate the d = 3 Ising model with the (finely-tuned) 
second and third neighbor interactions; see also the studies j^l 10,[ll| in the lattice- field- 
theory context. Their consideration could be viewed as an interesting application of the 



Monte Carlo renormalization group 



12j to exploit the virtue of the fixed point. (The Monte 



Carlo renormalization group provides an explicit realization of the renormalizat ion-group 
idea in the real space.) The aim of this paper is to develop an alternative approach to 
the elimination of corrections via the transfer-matrix method, and make the best use of its 
merits and characteristics. In fact, the four-spin interaction, appearing in our Hamiltonian 
(JIJ, is readily tractable with the transfer-matrix method, whereas it can make a conflict 
with the Monte Carlo simulation; in fact, the four-spin interaction disables the use of cluster 
update. (Probably, as for the Monte Carlo simulation, it might be more rewarding to 
enlarge the system size rather than incorporate extra interactions.) Second, the extended 
interactions appearing in our Hamiltonian (0) are taken from the proposal by Ma Q], who 
investigated the d = 2 Ising model and its renormalization-group flow. We consider that his 
renormalization-group scheme for d = 2 is still of use to our d = 3 well. Actually, in 

our transfer-matrix treatment, the system size along the transfer- matrix direction is infinite, 
and the remaining d = 2 fluctuations are responsible for the finite-size corrections. We 
demonstrate that Ma's scheme leads to satisfactory elimination of finite-size corrections in 
d= 3. 

The rest of this paper is organized as follows. In Sec. |Hj we explain the real-space 
decimation (renormalization group) for the d — 3 Ising model (0), and search for its fixed 
point. In Sec. IIII1 we perform extensive transfer-matrix simulation around this fixed point. 
Here, we utilized an improved transfer-matrix method, which is explicated in Appendix [XJ 
The last section is devoted to summary and discussions. 
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II. SEARCH FOR A SCALE-INVARIANT POINT: A REGIME OF ELIMINATED 
IRRELEVANT INTERACTIONS 



In this section, we search for a point of eliminated finite-size corrections of the extended 
Ising model, Eq. For that purpose, we set up a real-space renormalization group, and 
look for the scale-invariant (fixed) point; the result is given by Eq. (J2J). 

To begin with, we set up the real-space renormalization group. We consider a couple of 
rectangular clusters with the sizes 2x2 and 4x4; see Fig. These clusters are labeled 
by the symbols S and L, respectively. (Because we utilize the transfer-matrix method, the 
system sizes perpendicular to these rectangles are both infinite.) Decimating out the spin 
variables indicated by the symbol • of the L cluster, we obtain a reduced lattice structure 
identical to that of the S cluster. Our concern is to find a "scale invariance" condition with 
respect to this real-space renormalization group. 

Before going into the explicit formulation of the renormalization group (fixed-point anal- 
ysis), we explain briefly how we simulated the thermodynamics of these clusters. As men- 
tioned above, we employ the transfer-matrix method. The transfer-matrix elements for the 
5* cluster are given by the formula, 

%iM*> = «&) 1+36 (^& ^SJ; w£g w££) 1+ *> ( 2 ) 

where the component W s ^ s * denotes the local Boltzmann weight for a plaquette, Eq. (|A3|) . 
and the spin variables {Si} and {T;} (i — 1 ~ 4) denote the spin configurations for both sides 
of the transfer-matrix slice. The component (• • • ) 1+3b originates in the plaquette interactions 
perpendicular to the transfer-matrix direction, whereas the remaining part (■ ■ -) 1+6 comes 
from the longitudinal ones. The parameter b controls the boundary-interaction strength. 
Note that irrespective of b, the periodic-boundary condition is maintained; namely, all 2 x 2 
spins remain equivalent as b varies. Such a redundancy is intrinsic to the L = 2 system. 
Here, we consider this redundant parameter as a freely tunable one. (For example, a naive 
implementation of the periodic-boundary condition for a pair of spins may result in such an 
interaction as H = —JSiS 2 — JS 2 Si = — 2JS , 1 S , 2 . Apparently, such a duplicated interaction 
is problematic. Possibly, the interaction — (1 + b)JSiS2 with a certain moderate parameter 
b should be a favorable one. Significant point is that the periodic boundary condition is 
maintained with b varied.) We found that the choice b = 0.4 is reasonable because of the 
reasons mentioned afterward. 
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Similarly to the above, we constructed the transfer matrix for the L cluster as, 

1<*J<4 

with the 4x4 spin configurations {SV,} and {Tjj} under the periodic boundary condition. 
In this case (L = 4), we have no ambiguity as to the boundary interaction. 

Based on the above-mentioned simulation scheme, we calculate the location of the 
renormalizat ion-group fixed point. We impose the following "scale-invariance" conditions, 

(SiS 2 ) s = (S 1 S 2 )l (4) 
{SiSa)s = (5 , iS , 4 )l (5) 

(S1S2S3S4) s = (S1S2S3S4) l- (6) 

Here, the symbol (. . . )s(L) denotes the thermal average for the S (L) cluster, and the ar- 
rangement of spin variables {Si} is shown in Fig. ^ We solve the solution of the above 
equations numerically, and found that a non-trivial solution does exist at, 

(Jnn, Jnnn, h) = (0.1089828666435, 0.0445777727956, -0.0065117950492). (7) 

The last digits may be uncertain due to the numerical round-off errors. The result 
is to be compared with that of the preceding Monte Carlo study (Jnn, Jnnn, J^rd) 
(0.1109,0.03308,0.01402) 



j, where the authors incorporated the third- neighbor interaction 
J 3r d and omitted J\j instead. 

Let us mention a few comments. First, in the next section, we confirm that the fixed 
point is indeed a good approximant to the phase-transition point. This fact indicates that 
the above renormalizat ion-group analysis is sensible. Moreover, we calculated the fixed point 
Jnn = 0.2243904423106 for the conventional Ising model (Jnnn, Jri) = (0, 0). We again see 
that this transition point is in agreement with a critical point J^ N = 0.22165455(3) deter- 
mined with the Monte Carlo method [14]. (Hence, the choice of the boundary-interaction 
strength b = 0.4 is justified.) Second, we stress that the above renormalizat ion group is not 
intended to obtain (quantitatively reliable) critical point nor the critical indices. The aim 
of the above analysis is to truncate out the irrelevant interactions. The detailed analysis 
on criticality is made with the subsequent finite-size-scaling analysis. In other words, our 
numerical approach consists of two steps, and the remaining step is considered in the next 
section. 
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III. FINITE-SIZE SCALING ANALYSIS OF THE CRITICAL EXPONENTS v 

AND y h 



In the preceding Sec. |H1 we determined the position of the renormalization-group fixed 
point, Eq. (J7J). In this section, around the fixed point, we survey the criticality of the 
temperature- driven phase transition. Namely, hereafter, we dwell on the one-parameter 
subspace, 



which contains the renormalization-group fixed point at Jnn = Jnn- We anticipate that cor- 
rections to scaling (influence of irrelevant operators) should be suppressed in this parameter 
space. 

Throughout this section, we employ an improved version of the transfer-matrix method 
(Novotny's method) (To avoid confusion, we stress that in the above section, we used 
the conventional transfer-matrix method.) A benefit of Novotny's method is that we are 
able to treat arbitrary (integral) number of spins, N = 5, 6, . . . , 15, constituting a unit 
of the transfer-matrix slice; note that conventionally, the number of spins is restricted to 
N — 4, 9, 16, We explicate this simulation algorithm in Appendix 1X1 

A. Eliminated corrections to scaling 

In Fig. we plotted the scaled correlation length £/L for Jnn and a variety of system 
sizes N = 5, 6, . . . , 15. We evaluated the correlation length £ with use of the formula 
£ = 1/ ln(Ai/A2) with the dominant (sub- dominant) eigenvalue Ai (A2) of the transfer matrix. 
As explained in Appendix El the linear dimension L is simply given by, 



with the number of spins N; see Fig. |HJ 

From Fig. |2J we see a clear indication of criticality at Jnn ~ 0.11; note that the inter- 
section point of the curves indicates a critical point. Afterward, we compare this result to 
that of the conventional Ising model (Jnnn, Jo) = (0,0) to elucidate an improvement of 
the scaling behavior. Here, we want to draw reader's attention to the point that we treated 
various system sizes N = 5,6, ... ,15 with the aid of the Novotny method. Actually, in Fig. 






(9) 
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121 we notice that a variety of system sizes are available. Apparently, such an extension of 
available system sizes is significant in the subsequent detailed finite-size-scaling analyses. 

In Fig. El we presented the scaling plot (J NN - J* NN )L l/v -^/ L for 12 < iV < 15. with the 
scaling parameters J^ N = 0.11059 and v = 0.6245 determined in Figs. ElandEl respectively. 
We see that the data collapse into a scaling function satisfactorily; actually, we can hardly 
observe corrections to the finite-size scaling. 

In the above, we presented an evidence that the corrections-to-scaling behavior is im- 
proved in the parameter space, Eq. (jHJ). Lastly, as a comparison, we provide the data for 
the conventional Ising model; namely, we set (Jnnn, Jo) = (0,0) tentatively. In Fig. HI we 
plotted the scaled correlation length £/L for various Jnn- Apparently, the data suffer from 
insystematic finite-size corrections. The data scatter obscures the position of critical point. 
(Nevertheless, we should mention that the data imply Jt JM ~ 0.2, which does not contradict 
a recent Monte Carlo result J* NN = 0.22165455(3) 
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B. Phase-transition point J 



NN 

In the above, we obtain a rough estimate for the phase-transition point J^n- m this 
section, we determine the transition point more precisely. In Fig. we plotted the approx- 
imate transition point J^ rAr (L 1 ,L 2 ) for (2/(L\ + L 2 )) 2 . Here, the approximate transition 
point denotes the intersection point of the curves £/L (Fig. |2J) for a pair of system sizes 
(Li,L 2 ) (5 < Ni < N 2 < 15). Namely, the following equation, 

L ^( L l)\j NN =J* NN (LuL 2 ) = L ^ L ^\jNN=J* NN (Li,L 2 ) > ( 10 ) 

holds. In Fig. we notice that the data exhibit suppressed systematic finite-size deviation. 
Namely, the insystematic data scatter is more conspicuous than the systematic deviation. 
The least-squares fit to these data yields the transition point, 

J* NN = 0.11059(52), (11) 

in the thermodynamic limit L — ► oo. 



In order to check the reliability, we replaced the abscissa scale with 



where we used \ jv — 1.5868(3) and uo = 0.821(5) reported in Ref. 14j. (In the next section, 



(2/{L 1 +L 2 )) 



[15]. 



we make a consideration on the abscissa scale.) Thereby, we arrive at J^ N = 0.11062(43), 
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which is consistent with the above result. (The error margin may come from purely statistical 
one.) We confirm that the choice of the abscissa scale is not so influential. 

We notice that the transition point J NN (jllj) and the renormalizat ion-group fixed point 
Jnn are in good agreement with each other. This fact confirms that the renormalizat ion- 
group analysis in Sec. El is indeed sensible. As mentioned in Sec. El we do not require 
fine accuracy as to the convergence of J NN and Jnn- The aim of the renormalization- 
group analysis is to search for a regime of eliminated corrections rather than to obtain 
the precise location of the fixed point. Detailed analysis on criticality is performed in 
the subsequent finite-size-scaling analysis as demonstrated in the next section. (Actually, 
tuning the boundary- interaction parameter b (see Sec. EJ, we could attain better agreement 
between J NN and Jnn- However, such a refinement does not affect the subsequent finite- 
size-scaling analysis very much.) 



C. Critical exponents v and 

In Sec. IIII Al we presented an evidence of eliminated finite-size corrections. Encouraged 
by this result, in this section, we evaluate the critical exponents v and yn with use of the 
finite-size-scaling method (phenomenological renormalization group) [16]. 

In Fig. El we plotted the approximate correlation- length critical exponent, 



(T T \ 1 (T IT Wl f ggM/jO , gtt(jg)/jg) 

v(L 1 ,L 2 ) = HM/Lz)/ In — / — 

V c/Jnn OJnn 



' (12) 

Jnn—J^ n {L\,L2) 



for 2/{L\ + L2) with 5 < Aq < N2 < 15. With use of the least-squares fit to these data, we 
obtain the estimate, 

v = 0.6245(28), (13) 

in the thermodynamic limit. The data in Fig. El exhibit appreciable systematic finite-size 
corrections. More specifically, the systematic deviation (~ 2%) is almost comparable to 
the insystematic data scatter. (This fact indicates that we cannot fully truncate out the 
irrelevant interactions within the parameter space (Jnn, Jnnn, Jo)-) In this sense, the above 
(extrapolated) value, Eq. (fT2)l. may contain a systematic (biased) error. Afterward, we 
make a few considerations on the extrapolation scheme. (Because our work is methodology- 
oriented, we supply the least-squares- fit result as it is.) 
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In Fig. [7J we plotted the approximate exponent y^, 
VhiLxM) =ln ( ^ / ^ J /(21n(L 1 /L 2 )), (14) 

\ ' " .... j JNN=Jff N (Li,L2) 

for 2/(Li + L 2 ) with 5 < iVi < A^ 2 < 15. In order to incorporate the magnetic field H, 
we added the Zeeman term, —H'Y Ji Si 1 to the Hamiltonian (JTJ. Rather satisfactorily, the 
data yn(Li, L 2 ) exhibit suppressed systematic corrections; the systematic deviation is almost 
negligible compared to the insystematic one. The least-squares fit to these data yields the 
estimate, 

y h = 2.4709(73), (15) 

in the thermodynamic limit. 

Provided by the above estimates v and yh, we obtain the following critical indices through 
the scaling relations; 

a = 0.1265(84) (16) 
p = 0.3304(48) (17) 
7 = 1.213(11). (18) 

Let us provide comparative results with an alternative extrapolation scheme. We replaced 
the scale of abscissa in Figs. ED and |TJ with (2/(L x + L 2 )) UJ ; here, we set the exponent 
uj = 0.821(5) reported in Ref. [l4|. (As mentioned below, this scheme may overestimate 
the amount of systematic finite-size corrections.) Accepting this abscissa scale, we arrive at 
v = 0.6216(34) and yt = 2.4694(90). These values appear to be consistent with the above 
ones within the error margins, confirming that the extrapolation scheme is not so influential. 

We argue the underlying physics of the abscissa scale (extrapolation scheme) in detail. 
In principle, the exponent u governs the dominant (systematic) finite-size corrections. On 
the other hand, in the present simulation, we are trying to truncate out such systematic 
corrections. Hence, in our data analysis, the usage of the exponent u would be problematic. 
We consider that the systematic corrections should obey the scaling law like L~^ a ^ with a 
certain effective exponent u e ff > uj, at least, in the regime of 5 < iV < 15. Namely, we 
suspect that the abscissa scale with the exponent uj leads to an overestimation of systematic 
corrections. Actually, a recent Monte Carlo simulation reports the estimates v = 0.63020(12) 
and yh = 2.4816(1) |l4j ]. Here, we notice that their v indicates a non-negligible deviation, 
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whereas the value of yh is in good agreement with ours. This fact confirms the above 
observation that v{L\, L2) exhibits appreciable systematic corrections, and the extrapolated 
value may contain a biased error. Possibly, the adequate exponent u e jf would be even larger 
than the value w e ff = 1 utilized in Fig. |UJ Nevertheless, for the sake of simplicity, we do 
not pursue this issue further, and supply the least-squares-fit result as it is. 

Lastly, we mention a recent extensive exact-diagonalization result by Hamer jl7j |. who ob- 
tained v = 0.62854(79) and yh = 2.482(10). He investigated the quantum d = 2 transverse- 
field Ising model, relying on the belief that the quantum d = 2 Ising model should belong 
to the same universality class as the d = 3 Ising ferromagnet. The quantum version has an 
advantage such that the Hamiltonian elements are sparse (few non-zero elements), and one 
is able to treat a large cluster size 6x6. Comparing our data with his results, we notice 
that they are almost comparable with each other. Actually, the error margin of our yh is 
even smaller than his result, although we treated the d = 3 ferromagnet directly. 



IV. SUMMARY AND DISCUSSIONS 



So far, it has been considered that the transfer-matrix method would not be very useful to 
the problems in d > 3 because of its severe limitation as to the tractable system sizes. In this 
paper, we demonstrated that the corrections-to-scaling behavior of the d = 3 Ising model 
(JIJ is improved by adjusting the coupling constants to the values of the renormalizat ion- 
group fixed point, Eq. (J7J). Actually, corrections to scaling in Figs. El and 01 are eliminated 
significantly as compared to those in Fig. 0]for the conventional Ising model. Moreover, we 
succeeded in treating a variety of system sizes N = 5, 6, . . . , 15 with the aid of the Novotny 
method (Appendix [XJ; note that with the conventional approach, the available system sizes 
are restricted to N = 4, 9, 16, ... . Apparently, such an extension of available system sizes 
provides valuable information on criticality. Owing to these improvements, we analyzed 
the criticality of the d = 3 Ising model with the transfer-matrix method, and obtained the 
critical indices v = 0.6245(28) and y h = 2.4709(73). 

As mentioned in Introduction, an attempt to eliminate finite-size corrections has been 

n n 

pursued |8j in the context of the Monte Carlo renormalizat ion group |12j. We consider 
that an approach with the transfer-matrix method is also of use because of the following 
reasons. First, we accepted a simple renormalization-group scheme shown in Fig. [0 As 
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mentioned in Introduction, this scheme was introduced originally as for the d = 2 Ising 
model Q]. The advantage of the transfer-matrix method is that the system size along the 
transfer-matrix direction is infinite, and the remaining d = 2 fluctuations are responsible for 
the finite-size corrections. Hence, such a (d — 2)-like renormalization group is still of use 
to achieve elimination of corrections satisfactorily. Second, the transfer-matrix method is 
capable of the four-spin interaction appearing in our Hamiltonian (0). On the other hand, 
the Monte Carlo sampling conflicts with such a multi-spin interaction, because the multi- 
spin interaction disables the use of cluster- up date algorithm. (Probably, an effort toward 
enlarging the system size would be rewarding from a technical viewpoint.) 

In addition to these merits, we would like to emphasize again the point that the transfer- 
matrix approach with the Novotny method allows us to treat a variety of system sizes 
N = 5, 6, . . . , 15. We consider that Novotny's method combined with the elimination of 
finite-size corrections would be promising to resolve the (seemingly intrinsic) drawback of 
the transfer-matrix method in d > 3. As a matter of fact, the basic idea of the present 
scheme would be generic, and it might have a potential applicability to a wide class of 
systems. An effort toward this direction is in progress, and it will be addressed in future 
study. 
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APPENDIX A: NOVOTNY'S TRANSFER-MATRIX METHOD 

We explain the details of the transfer-matrix method utilized in Sec. IIHI (To avoid 
confusion, we remind the reader that in Sec. [HI we utilized the conventional transfer-matrix 







which enables us to 



method.) Our method is based on Novotny's formalism 
consider an arbitrary (integral) number of spins v A r , constituting a unit of the transfer- 
matrix slice even for d > 3; note that conventionally, the number of spins is restricted to 
AT = 4, 9, 16, ... . We made a modification to the Novotny formalism in order to incorporate 
the plaquette-type interactions. We already reported this method in Ref. p, where we 
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studied the multicriticality of the extended d = 3 Ising model f^Lsj] . In the present paper, 
we implemented yet further modifications such as Eqs. (jA8|) - (jA10|) . Hence, for the sake of 
self-consistency, we explicate the full details of the simulation scheme. 

Before going into details, we mention the basic idea of the Novotny method. In Fig. 
IH1 We presented a schematic drawing of a unit of the transfer-matrix slice. Note that in 
general, a transfer-matrix unit for a <i-dimensional system should have a (d— l)-dimensional 
structure, because it is a crosssection of the ci-dimensional manifold. However, as shown in 
Fig. |Sl the constituent N spins form a d = 1 (coiled) alignment rather than d = 2. The 
dimensionality is raised effectively to d = 2 by the y/~N th-neighbor interactions among these 
N spins; This is the essential idea of the Novotny method to constitute a transfer-matrix 
unit with arbitrary number of spins even for d = 3. 

In the following, we present the explicit formulas for the transfer-matrix elements. We 
decompose the transfer matrix into the following three components: 

T(v) = T (leg) T iplanar \v) T {run9 \v), (Al) 

where the symbol denotes the Hadamard (element by element) matrix multiplication. 
Note that the product of local Boltzmann weight gives rise to the global one. The physical 
content of each component is shown in Fig. |H| 

The explicit expression for the element of T^ le9 ^ is given by the formula, 

T (leg) _ IftAUK _ W S(j,l)S(j,2) w S(j,2)S(j,3) w S(j,N)S(j,l) , A „v 

1 ij - W A \J) - VV S(i,l)S(i,2) W S(i,2)S(i,3) ■ ■ ■ W S{i,N)S{i,l) ' 

where the indices i and j specify the spin configurations of both sides of the transfer-matrix 
slice. More specifically, the spin configuration {S(i, 1), S(i, 2), . . . , S(i, N)} is arranged along 
the leg; see Fig. |HJ The factor W s ^ s 4 denotes the local Boltzmann weight for the plaquette 
spins {Si} (i — 1 ~ 4); 

Wff = exp (- f-^pCS'iSa + S 2 S, + S,S 3 + SM - ^^{S l S 4 + S 2 S 3 ) - ^S 1 S 2 S 3 S 4 

(A3) 

Notably enough, the component T^ 9 ^ is nothing but a transfer matrix for the d = 2 Ising 
model. The remaining components T^ lanar ^ and T^ run9 ^ introduce the v^th -neighbor cou- 
plings, and raise the dimensionality effectively to d — 3. 
The component T^ planar ^ is given by, 

T \f anar \v) = (i\AP v \t), (A4) 
12 



with 

v = y/N, (A5) 

where the matrix P denotes the translation operator. Namely, the state P\i) represents 
a shifted configuration {S(i,m + 1)}. Hence, the insertion of P^ introduces the x/iVth- 
neighbor interactions among the N spins [lj . Similarly, we propose the following expression 
for the component T^ run9 ^; 



j. In the following, we propose 



T irun 9){v) = ^ ^ B {{py)) ^ j (A6) 

with, 

ail ® oi) b (I*) ® io) = fi ^ a? ) 

m=l 

The meaning of the formula would be apparent from Fig. 

The above formulations are already reported in Ref. 
a number of additional improvements. First, we symmetrize the transfer matrix with the 
replacement 

T{v) -^T{v) QT(-v). (A8) 

Correspondingly, we substitute the strength of the coupling constants J a — > J a /2 in order 
to compensate the above duplication. Apparently, with the symmetrization, the symmetry 
of descending (m = N, N — 1, . . .) and ascending (m = 1,2,...) directions is restored. 
Moreover, we implement the following symmetrizations, 

(i\AP v \i) -> (i\AP v \i)(i\p- v A\i), (A9) 

and, 

((i\®(j\)B((py))®(p v \j)))^ 

® B ((P*\i)) 9 (P v \j))) {((i\P- v ) ® ((j\P- v )) B(\i)® , (A10) 

as to Eqs. (|A4|) and (|A6|1 . respectively. Again, we have to redefine the coupling constants 
to compensate the duplication. 
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FIG. 1: A schematic drawing of our real-space renormalization group (decimation) for the d = 3 
Ising model with the extended interactions, Eq. (|T]). As indicated, the thermodynamics is simulated 
with the transfer-matrix method. Imposing a criterion of scale invariance, Eq. @, we determined 
the renormalization-group fixed point, Eq. ( [7| ) : numerically; see text for details. 
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FIG. 2: Scaled correlation length £/X is plotted for the nearest-neighbor interaction Jnn and N = 
5, 6, . . . , 15 (N = L 2 ); note that we survey the parameter space (JHJ) including the renormalization- 
group fixed point (J7J). We observe a clear indication of criticality at J^tv ~ 0.11. Apparently, the 
finite-size-scaling behavior is improved as compared to that of Fig. |1] for the conventional Ising 
model. 
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FIG. 3: The scaling plot for the correlation length, (Jnn ~ Jnn)^ 1 ^-^/!^, is shown for the system 
sizes (+) N = 12, (x) 13, (*) 14, and (□) 15; note the relation N = L 2 . Here, we accepted the 
scaling parameters, = 0.11059 and v = 0.6245, determined in Figs. ElandEl respectively. We 
again confirm that corrections to scaling are suppressed significantly. 





FIG. 4: Tentatively, we turned off the extended interactions (Jnnn = and J\j = 0), and 
calculated the scaled correlation length £/L for various Jnn and iV = 5, 6, . . . , 15. We notice that 
the data are scattered as compared to those in Fig. |2j 



16 



0.114 
0.113 - 
0.112 - 
0.111 - 
0.11 - 
0.109 
0.108 h 
0.107 







+ 



+ + 



+ ++ + 
+ 



+ 



0.05 



0.1 

(2/(L 1 + L,)) : 



0.15 



0.2 



FIG. 5: The approximate critical interaction J^ N (Li, L 2 ) is plotted for (2/(Li + L 2 )) 2 with 
5 < Ni < N 2 < 15 (Li ;2 = ^N 1>2 ). The least-squares fit to these data yields J^ N = 0.11059(52) 
in the thermodynamic limit L — > 00. 
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FIG. 6: The approximate correlation-length critical exponent v(Li,L 2 ) is plotted for 1j{L\ + L 2 ) 
with 5 < iVi < N 2 < 15 (L 1<2 = ^/W^)- The least-squares fit to these data yields v = 0.6245(28) 
in the thermodynamic limit L —> 00. 
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FIG. 7: The approximate critical exponent yf l {L\,L2j is plotted for 2/{L\ + L-z) with 5 < N\ < 
N% < 15 (Li t 2 = y/Nip)- The least-squares fit to these data yields y^ = 2.4709(73) in the 
thermodynamic limit L — > 00. 




FIG. 8: A schematic drawing of a unit of the transfer-matrix slice for the d = 3 Ising model with 
the extended interactions (JTj). The contributions from the "leg," "planar," and "rung" interactions 
are considered separately; see Eq. (|Alj) . Within the transfer-matrix slice, the arrangement of the 
constituent spins is one-dimensional (coiled structure). The dimensionality is raised up to d = 2 
by the bridges between the y/Nth- neighbor spins along the leg. 
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